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We discuss the extraction of the electromagnetic and spin polarisabilities of nucleons from lattice 
QCD. We show that the external field method can be used to measure all the electromagnetic and 
spin polarisabilities including those of charged particles. We then turn to the extrapolations required 
to connect such calculations to experiment in the context of finite volume chiral perturbation theory. 
We derive results relevant for lattice simulations of QCD, partially-quenched QCD and quenched 
QCD. Our results for the polarisabilities show a strong dependence on the lattice volume and quark 
masses, typically differing from the infinite volume limit by ~ 10% for current lattice volumes and 
quark masses. 
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Compton scattering at low energies is an invaluable tool with which to study the electromagnetic structure of 
hadrons. At very low photon energies, the Compton amplitude is dominated by point-like photon scattering from 
the total charge and magnetic moment of the target hadrons. As the frequency increases, contributions beyond 
point-like scattering enter and one begins to resolve the hadronic response to an applied electromagnetic field. For 
unpolarised scattering on spin one-half objects, the first structure dependent contributions in this energy expansion of 
the amplitude are the electric polarisability, a, and the magnetic polarisability, (3. These quantities reflect the ability 
fvq , of the hadron's components to align or anti-align themselves in response to an applied electric or magnetic field. For 
the proton and neutron, the positivity of the accepted experimental values of these polarisabilities (a p — 12.0 ± 0.6, 
(3 P = 1.9 T 0.6, On = 12.5 ± 1.7 and [3 n = 2.7 T 1.8 in units of 10~ 4 fm 3 Q) indicates that both nucleons are 
diamagnetic objects. Recent experimental advances 0,0 have also allowed the extraction of certain combinations of 
target polarisation-dependent observables in Compton scattering. These involve the so-called spin polarisabilities , 
conventionally labeled 71-74, and they have consequently been investigated in numerous theoretical and further 
experimental studies. Although the classical interpretation of spin-dependent Compton scattering is less clear, the spin 
, polarisabilities encode additional fundamental properties of the nucleon. Compton scattering observables, however, 
are not limited to these six parameters. Higher order quasi-static properties of the nucleon appear from further terms 
in the energy expansion of the amplitude. These higher-order polarisabilities [J], as well as generalised polarisabilities 
Si ■ @ (which arise in the singly (doubly) virtual Compton scattering process, 7* A — > 7*-*'A) allow for an even finer 
resolution of the electromagnetic structure of hadrons at low energies. 

While experimentally one is hoping to open further windows through which to view hadronic electromagnetic 
structure, theoretically one ultimately hopes to understand how hadronic polarisabilities arise from the basic elec- 
tromagnetic interaction of the photon with quarks that are bound to form the hadrons. The electric and magnetic 
polarisabilities should naively scale with the volume of the hadron. However, this expectation overestimates the ob- 
served polarisabilities by four orders of magnitude, indicating that the nucleon's constituents are strongly coupled. 
Lattice techniques provide a method to investigate the non-perturbative structure of hadrons directly from QCD. 
In particular, the various hadron polarisabilities can be computed. Comparison of these results with experimental 
determinations would provide stringent tests of the lattice method's ability to reproduce the structure of physical 
hadronic states; for the individual spin polarisabilities that have not been measured, the lattice approach may be the 
only way to determine them. On the lattice, direct calculations of the required hadronic current-current correlators 
are difficult and so far have not been attempted. However significant progress has been made 0, 0, in extracting 
the electric and magnetic polarisabilities by performing quenched lattice calculations in constant background electric 
and magnetic fields respectively and studying the quadratic shift in the hadron mass that is induced (essentially an 
application of the Feynman-Hellman theorem) . These studies have investigated the electric polarisabilities of various 
neutral hadrons (in particular, the uncharged vector mesons and uncharged octet and decuplet baryons), and the 
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magnetic polarisabilities of the baryon octet and decuplet, as well as those of the non-singlet pseudo-scalar and vector 
mesons. As we shall discuss below, generalisations of these methods using non-constant fields allow the extraction 
of the spin polarisabilities from spin-dependent correlation functions and also allow the electric polarisabilities to 
be determined for charged hadrons. More generally, higher-order polarisabilities and generalised polarisabilities are 
accessible using this technique. 

As with all current lattice results, these calculations have a number of limitations and so are not physical predictions 
that can be directly compared to experiment. For the foreseeable future, lattice QCD calculations will necessarily use 
quark masses that are larger than those in nature because of limitations in the available algorithms and computational 
power. Additionally, the volumes and lattice-spacings used in these calculations will always be finite and non- vanishing, 
respectively. For sufficiently small masses and large volumes, the effects of these approximations can be investigated 
systematically using the effective field theory of the low energy dynamics of QCD, chiral perturbation theory (xPT) 
lift fill]. 1 In this paper we shall perform an analysis of the nucleon electromagnetic and spin polarisabilities at 
next-to-leading order (NLO) in the chiral expansion. We do so to discuss the infrared effects of the quark masses 
and finite volume in two-flavour QCD and its quenched and partially-quenched analogues (QQCD and PQQCD). The 
polarisabilities of the hadrons are particularly interesting in this regard since they are very sensitive to infrared physics 
and their quark mass and finite volume dependence is considerably stronger than that expected for hadron masses 
and magnetic moments. This should be physically evident given that the polarisabilities scale with the volume. In 
essence, chiral perturbation theory provides a model independent analysis of the modification of the nucleon's pion 
cloud in a finite volume. When the charged pion cloud is influenced by to the periodic boundary conditions imposed 
on the lattice, the nucleon's response to external electromagnetic fields is altered compared to that at infinite volume, 
and in most cases the effects are dramatic. A particularly striking oddity that we find in this analysis is a modification 
of the Thompson cross section at finite volume. This can be explained through the physics of chiral loop corrections 
to point-like hadron structure. 

If future lattice QCD simulations are to provide physical predictions for the electromagnetic and spin polarisabilities, 
careful attention must be paid to both the chiral and infinite volume extrapolations. To illustrate this point, we present 
our results at representative values of the quark mass, finding significant effects. We also use our quenched chiral 
perturbation theory results to assess the volume dependence of the quenched data at the lightest pion masses used 
in Refs. 0,13- While the quenched theory contains unphysical low energy constants (LECs) and the convergence of 
the chiral expansion is questionable at these pion masses, we can still provide an estimate of the volume dependence 
of quenched data for the nucleon polarisabilities using our results. Such an estimate is achievable because the 
corresponding polarisabilities in the unquenched theory do not depend on phenomenologically undetermined LECs 
at the order of the chiral expansion to which we work. At the lightest quark masses used in the existing quenched 
lattice simulations, ~ 0.5 GeV, we find strong sensitivity to the lattice volume (as large as 10%). The effects will 
only increase as the pion mass is brought closer to that in nature. Clearly careful chiral and volume extrapolations of 
polarisabilities are mandated to connect lattice calculations to real world QCD. 

To begin our investigation of nucleon polarisabilities in lattice QCD, we discuss in Sec.|n]the kinematics of Comp- 
ton scattering and define the electromagnetic and spin polarisabilities that are the primary focus of this work. In 
Sec. IIIII we perform a general analysis of the external field method pertaining to all electromagnetic and spin polar- 
isabilities. We discuss how suitable background fields can be used in lattice QCD simulations to determine the spin 
polarisabilities and, more generally, generalised polarisabilities (though we limit our discussion of these in the present 
paper). Following this we introduce the low-energy effective theories of QCD (%PT), quenched QCD (QxPT) and 
partially-quenched QCD (PQxPT). These effective theories provide the model independent input necessary for cal- 
culating the quark mass and lattice volume dependence of polarisabilities. We focus primarily on PQxPT in Scc. lIVI 
discussing the relation to X-PT where relevant, and relegating the peculiarities of QxPT to Appendix IA1 Our results 
for the dependence of the nucleon polarisabilities on quark masses and the lattice volume are presented in Sec.0 We 
provide detailed plots relevant for full QCD simulations of polarisabilities showing the dependence on quark masses 
and lattice volumes. We also estimate the quenched QCD volume dependences of the polarisabilities at a pion mass 
typical of existing quenched lattice data. Very small volumes are discussed in Appendix ]B\ and a glossary of finite 
volume functions required to evaluate the polarisabilities in a periodic box appears in Appendix ICl Lastly, Sec. I VII 
consists of a concluding discussion of our results. 



The effects of the lattice discretisat ion are short d istance in nature, and while some of them can be analysed in an extension of the 
effective field theory described here 1 1 d. 1 1 1 ll~3 . 1 1 r J , others are not accounted for. Here we will assume that a continuum extrapolation 
has been performed. 
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II. COMPTON SCATTERING ON SPIN-HALF HADRONS 



The real Compton scattering amplitude describing the elastic scattering of a photon on a spin-half target such as 
the proton or neutron can be parameterised as 



T 7JV = Ai(w,0) e' ■ e + A 2 (uj 7 9) e'-fee - k' + i A 3 (u, 6) a ■ (e" x e) + i A 4 (uj 7 9) a ■ (k' x k) e' ■ e 

+iA 5 (u,6)a- \(e' x k)e- k> - (ex k')e' ■ k] +iA 6 (u,9)a- [(?' x k')e- k> - (ex k) e ' ■ k 



(1) 



where we have chosen to work in the Breit frame of the system and the incoming and outgoing photons have momenta 
k = (uj, k = oj k) and k' = (uj, k! = u> k'), and polarisation vectors e and e', respectively. The Ai(u>, 9), i = 1 ... 6, are 
scalar functions of the photon energy and scattering angle, cos9 = k-k . It is convenient to work in Coulomb gauge 
throughout where eo = e' = (the physical amplitudes are gauge invariant). 

The functions, Ai, determining the Compton amplitude can be separated into a number of pieces. The Born terms 
describe the interaction of the photon with a point-like target with mass, Mjv, charge, e Z (where e > 0), and magnetic 
moment, /i. These terms reproduce the Thomson-limit an d q uadratic frequency pieces |l6l of unpolarised scattering 
and the Low-Gell-Mann-Goldberger low energy theorems jl7t llSl ] for spin-dependent scattering. The remaining parts 
of the amplitude describe the structural response of the target. Expanding the amplitude for small photon energies 
relative to the target mass and keeping terms to 0(ui 3 ) one can write 
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describing the target structure in terms of the electric, magnetic and four spin polarisabilities, a, /3, and 71-4, 
respectively. In the conventions above, the spin polarisabilities receive contributions from the anomalous decay 
7T° — > 77 (shown in Fig.|21below). This contribution varies rapidly with energy and is omitted from the polarisabilities 
in some conventions. Higher order terms in the energy expansion can be parameterised in terms of higher-order 
polarisabilities Q . The more general process of virtual (and doubly- virtual) Compton scattering at low energies can 
similarly be described in terms of generalised polarisabilities 0. We will focus in the six polarisabilities defined above. 

The goal of this paper is to determine the quark mass and volume dependence of the polarisabilities defined above 
to allow accurate extraction of their physical values from lattice calculations. Before we do this we shall discuss how 
these lattice calculations may be implemented. 



III. COMPTON SCATTERING AND POLARISABILITIES ON THE LATTICE 



Lattice QCD provides a way to study the polarisabilities of hadrons from first principles. There are two ways to do 
this. The method most reminiscent of the experimental situation is to study the (Euclidean space) four point Green 
function defining the Compton scattering tensor directly (the photon fields are amputated). By measuring the large 
Euclidean time behaviour of this correlator, the hadron matrix elements of the two vector currents can be extracted. 
In principle, by calculating particular Lorentz components of the Compton tensor with various different source and 
sink spin states, all six electromagnetic and spin polarisabilities and their higher order and generalised analogues 
can be extracted. However, this is a complicated task, requiring the evaluation a large number of quark propagator 
contractions resulting from quark-line disconnected diagrams which are statistically difficult to determine. At present 
this approach is too demanding for the available computational resources and has not been attempted. 

The second method is based on measuring the response of hadronic states to fixed external fields. A number of 
exploratory quenched QCD studies have been performed in this approach. The pioneering calculations of Refs. [T9I 
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l2Ci l2lU22U 23ll24| attempted to measure the nucleon axial couplings, magnetic moments and electric dipole moments 
by measuring the linear shift in the hadron energy as a function of an applied external weak or electromagnetic field. 
As discussed in the Introduction, various groups |E S H3 have also used this approach to extract electric and 
magnetic polarisabilities in quenched QCD by measuring a quadratic shift in the hadron energy in external electric 
and magnetic fields. The method is not limited to electroweak external fields and can be used to extract many matrix 
elements such as those that determine the moments of parton distributions and the total quark contribution to the 
spin of the proton |26| . Here we focus on the electromagnetic case. 

The Euclidean space (x^ = r) effective action describing the gauge and parity invariant interactions of a non- 
relativistic spin-half hadron of mass M and charge q with a classical U(l) gauge field, A^(x, r), is 

S cS [A] = J d 3 xdr C eS (x, r; A) , (3) 

for the Lagrangian 



C e s(x, r; A) = &(x,t) 



2m (-iE lEl a -ExE + y MlMl ff -HxH + lMlE2 a i E l] H 3 + lEl M 2 a t H t3 E^ 



tf(2,r) + ... , (4) 



where E — —4^A(x,t) — VA±(x, t) and H = V x A(x, r) are the corresponding electric and magnetic fields, X — 
-§^X denotes the Euclidean time derivative, — ^{d l X^ + d^X 1 ), and the ellipsis denotes terms involving higher 
dimensional operators. By calculating one- and two-photon processes with this effective Lagrangian, it is clear that 
the constants that appear in Eq. Q are indeed the relevant magnetic moment and electromagnetic and multipolc 
polarisabilities [23] [these are simply related to the polarisabilities defined in the previous section as: 7s 1 e 1 — —(71 + 
73), 'Jm 1 m 1 — 74, 7_Ei m 2 — 73 and Jm 1 e 2 = 72+74]- The Schrodinger equation correpsonding to Eq. 0} determines the 
energy of the particle in an external U(l) field in terms of the charge, magnetic moment, and polarisabilities. Higher 
order terms in Eq. Q (which contain in part the higher order polarisabilities Q) can be neglected for sufficiently 
weak external fields. For a magnetic field, the minimally coupled terms generate towers of Landau levels and for a 
constant electric field the same terms accelerate charged particles. 

Lattice calculations of the energy of a hadron in an external U(l) field are straight-forward. One measures the 
behaviour of the usual two-point correlator on an ensemble of gauge configurations generated in the presence of the 
external field. This changes the Boltzmann weight used in selecting the field configurations from det [Ip+ m] exp [— S g ] 

to det lp+ iQA + m exp [— S g ], where If) is the SU(3) gauge covariant derivative, Q is the quark electromagnetic 

charge operator, and S g is the usual SU(3) gauge action. Since calculations are required at a number of different values 
of the field strength in order to correctly identify shifts in energy from the external field, this is a relatively demanding 
computational task (although it is at least conceptually simpler than studying the four-point function). In general 
one must worry about the positivity of the fermionic determinant calculated in the presence of a back grou nd field, 
however for weak fields, positivity is preserved. The exploratory studies of Refs. |rj ItI IsL liH I20L l2ll I22I I23LI24I I25I ] used 
quenched QCD in which the gluon configurations do not feel the presence of the U(l) field as the quark determinant 
is absent. In this case, the external field can be applied after the gauge configurations had been generated and is 
simply implemented by multiplying the SU(3) gauge links of each configuration by link variables corresponding to the 
fixed external field: {U£(x)} — > {U£(x) exp[i e a A^]}, where a is the lattice spacing. These studies are interesting 
in that they provide a proof of the method, however the values of the polarisabilities extracted have no connection to 
those measured in experiment. 

It is clear from Eq. Q that all six polarisabilities can be extracted using suitable space and time varying background 
fields if the shift of the hadron energy at second order in the strength of the field can be determined. One can also see 
this because the Compton tensor appears explicitly as the second-order connected term in the expansion of hadronic 
two-point correlation function in weak background fields [2(J . Previous studies 0> H> US have employed constant 
electric and magnetic fields to determine the corresponding polarisabilities in quenched QCD. Here we perform a 
more general analysis to show how the spin polarisabilities and the electric polarisabilities of charged particles can be 
obtained. 

In order to determine the polarisabilities, we consider lattice calculations of the two-point correlation function 

C ss ,(p,t;A)= [ d 3 xe^ s (0\ Xs (x,T)xl(0,0)\0)A, (5) 
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where Xs (x, r) is an interpolating field with the quantum numbers of the hadron under consideration (we will focus 
on the nucleons) with z component of spin, s, and the correlator is evaluated on the ensemble of gauge configurations 
generated with the external field, . 

For uncharged hadrons at rest in constant electric and magnetic fields, it is simple to show that this correlator falls 
off exponentially at large times with an energy given by the appropriate terms in Eq. (@J owing to the constancy of 
the effective Hamiltonian. However for space-time varying fields, charged particles or states of non-zero p, a more 
general analysis is needed. This is most easily formulated using the effective field theory (EFT) defined by Eq. 
For weak external fields (such that higher order terms in Eq. (0J can be safely neglected) , the small p and large r 
dependence of this QCD correlation function is reproduced by the equivalent correlator calculated in the effective 
theory corresponding to the Lagrangian, Eq. Q . That is 



C ss ,(p,t;A) = 



d 3 xe ii! - s 



1 



Z oS [A] 



V&TW r)* T s , (0, 0) exp (S eS [A\) , 



(6) 



where -Z e ff[A| = J "D^D^f exp (— 5 c ff[A]). Since the right-hand side of Eq. © is completely determined in terms of 
the charge, magnetic moment and polarisabilities that we seek to extract, fitting lattice calculations of C ss i (p, r; ^4) in 
a given external field to the effective field theory expression will enable us to determine the appropriate polarisabilities. 
In the above equation we have assumed that the ground state hadron dominates the correlator at the relevant times. 
For weak fields this will be the case. However one can consider additional terms in the effective Lagrangian that 
describe the low excitations of the hadron spectrum that have the same quantum numbers as the hadron under study. 
This will lead to additive terms in Eq. 10 that depend on the mass, magnetic moment and polarisabilities of the 
excited hadron instead of those of the ground state. With precise lattice data, the properties of these excited states 
can also be determined. 

In many simple cases such as constant or plane-wave external fields, the EFT version of C ss i(p, t; A) can be 
determined analytically in the infinite volume, continuum limit [2^. However in finite lattice spacing and at finite 
volume, calculating C ss '(p, t;A) in the EFT becomes more complicated. In order to determine the EFT correlator, 
we must invert the matrix K. defined by 



S\ at t[A] = ^^^^l(x,T x )K, ss '[x,T x ,y,T y ;A]^ s (y,T y ), 



(7) 



x,r x y,r y s.s' 



where <Si a tt[^4] is a discretisation of the EFT action in which derivatives are replaced by finite differences (the time 
derivative is given by a forward difference as we can ignore anti-particles). K. has dimension 47V 2 where TV; is the 
number of lattice sites. For the most general space-time varying external field, this must be inverted numerically; 
given a set of lattice results for the correlator, Eq. © is repeatedly evaluated for varying values of the polarisabilities 
until a good description of the lattice data is obtained. 

For weak fields such that \A^(x, r)| 2 <g; Aq CD for all x and t, a perturbative expansion of /C _1 in powers of the 
field can be used. This corresponds to the series of diagrams in Fig. To extract all six polarisabilities using such 
an analysis, we need to consider a number of different fields; lattice calculations of the correlators in Eq. © using 
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for a number of different choices for the strength parameters, a, and hi (with \di\, <C ^qcd)' are sufficient to 
determine the full set of polarisabilities. 2 By measuring correlators for different spin configurations (including those 
that flip spin), we can reduce the number of fields required to extract the polarisabilities. 

As an example, the behaviour of the correlator in the field (x) (which corresponds to a constant electric field 
in the x\ direction) is given by 



C ss ,(p,t;A {1) ) = 5 SiS ,exp{-^ [a x (q 2 r 2 + UMna) - 3iqr Pl ] } e - M - e -« 2 + {a\) 



exp 



-(M + 2naa\)T 



1 2 

6M 



+ 0(al). 



(9) 



In this case, the perturbative series has been resummed exactly in the continuum, infinite volume limit and the higher 
order corrections come from terms omitted in Eq. Q). For electrically neutral particles, the exponential fall-off of this 
correlator determines the polarisability a once the mass M has been measured in the zero-field case. When a charged 
particle is placed in such a field it undergoes continuous acceleration in the x\ direction (this is described by the r 3 

term in the exponent). However at times small compared to "^ 6 Jf , the correlator essentially falls off exponentially. 
Matching the behaviour of Eq. © to lattice data for a charged hadron will again enable us to determine the electric 
polarisability, a. 

As a second analytic example, we consider one of the multipole polarisabilities. In the presence of the field A^(x), 
which corresponds to a more complicated electric field E^(x) — (— r, z^-,0), we find that 



C^(p,t]A (6) ) 
Cu(p,t;A {6) ) 



exp 



2tt 



■ a 6 b 6 j El E! t 



(10) 



independent of p and the ellipsis denotes terms cubic in the field that have been neglected in Eq. (0J . Whilst the 
individual correlators, Cjj and Cj_j_ , have relatively complicated time-momentum behaviour involving q and a as well 
as 7_e 1 b 1 , this becomes very simple in the ratio and "/e 1 e 1 can be determined cleanly. 

Analogous results can be derived for the other fields in Eq. ijHJ), however to take into account the finite lattice spacing 
and periodic finite volume nature of the underlying lattice simulations to which the EFT description is matched, the 
correlator is most easily calculated by inverting the matrix K. numerically. This also allows for more general choices of 
fields. If we seek to extract higher order polarisabilities, the Lagrangian in Eq. must be extended to include higher 
dimension operators 0. At this order, relativistic corrections and three-photon couplings also need to be included. 
Correlation functions similar to those in Eq. JSJ involving two different external momenta will allow us to also extract 
the generalised polarisabilities 



IV. HEAVY BARYON X PT 



To calculate the quark mass and volume dependence of the nucleon polarisabilities, we use heavy baryon chiral 
perturbation theory (HB^PT) as was first constructed in Refs. (29l. l30t l3lL l3^ | . In current lattice calculations, valence 
and sea quarks are often treated differently, with sea quarks either absent (quenched QCD) or having different masses 
than the valence quarks (partially-quenched QCD). 3 The extensions of HB^PT to quenched HB%PT [3^, |3(| and 
partially quenched HB%PT 37, 38] to accommodate these modifications are also well established and have been used 
to calculate many baryon properties. In this section, we will primarily focus on the two flavour partially-quenched 
theory and briefly introduce the relevant details following the conventions set out in Ref. j^. Since QCD is a special 
limit of the partially-quenched theory, our discussion also encompasses two flavour xPT- Additional complications in 
quenched xPT are relegated to Appendix lAl 



2 These fields correspond to real E and H fields in Minkowski-space for real- valued a; and bi. Since pe riodic s patial boundary conditions 
are envisaged for the link variables, there are quantisation conditions that must be satisfied by the a; |2(1|2I| . For example, qi<i2 = 2 ™£ 
for each of the quark charges qi. The more complicated fields in Eq. require two parameters to satisfy these conditions. 

3 At finite lattice spacing, different actio ns can even be used for the different quark sectors (e.g., staggered sea quarks and domain wall 
valence quarks). As was shown in Refs. Il "ILi lF3 . the lattice spacing corrections to baryon electromagnetic properties are expected to 
be small, as they can not enter at tree level, and for current simulations with o,Aq CD ~ m q , they generally enter at leading loop order 
through valence-sea meson masses. In our work we assume a continuum extrapolation has been performed. 
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Pseudo-Goldstone mesons 



We consider a partially-quenched theory of valence (u, d), sea (J, I) and ghost (u, d) quarks with masses contained 
in the matrix 



m Q = diag(m u , m^m^mi, m fi ,mj) 



(11) 



where 



*u,d 



m Uj d such that the path-integral determinants arising from the valence and ghost quark sectors exactly 



cancel. The corresponding low-energy meson dynamics are described by the PQxPT Lagrangian. At leading order 
this is given by 



= £- s tr [V^Vf.Y,] + \—str 
8 ' 4 



where the pseudo-Goldstone mesons are embedded non-linear ly in 



with the matrix $ given by 
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and 



M = 
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exp 



M x r 
X M 



f 



(12) 



(13) 



(14) 



J du 



"dd 



'dl 



(15) 



and where 4>o = str ($) /v2- The upper left 2x2 block of M corresponds to the usual valence-valence mesons, the 
lower right to sea-sea mesons and the remaining entries of M to valence-sea mesons. Mesons in M are composed of 
ghost quarks and ghost anti-quarks and thus bosonic. Mesons in \ contain ghost-valence or ghost-sea quark-anti- 
quark pairs and are fermionic. In terms of the quark masses, the tree-level meson masses are given by 



A 



( m Q)u + ( m Q)i 



(16) 



where Q = (it, d, j, I, u, d). The terms proportional to a$ and mo in Eq. 112|) involve the flavour singlet field and are 
only relevant in the quenched theory (see Appendix 0; in both PQxPT and xPT the singlet meson acquires a large 
mass through the strong U(l)^i anomaly and can be integrated out, leading to a modified flavour neutral propagator 
that contains both single and double pole structures [3!j. 

In the above Lagrangian, we have minimally coupled electromagnetism (the U(l) gauge field is again denoted by 
and its field strength tensor F^ v — d^A u — d v A^) to the theory through the chiral, and U(l) gauge covariant 
derivative 



with the vector current 



V 



[ £ (0" - ieeA")^ + tHW-ieQA^S] 



(17) 



(18) 



depending on the quark charge matrix, Q. In coupling electromagnetism to this theory, we must specify how the 
quark charges are extended to the partially-quenched theory. We choose: 



di&g(q u ,q d ,qj,qi,q u ,q d ) , 



(19) 



though other arrangements are possible. However, one must set qj +qi ^ in order to retain sensitivity to the full set 
of LECs that appear in two flavour xPT E3,EjJ. In addition to the Lagrangian, Eq. I|12[l. the anomalous couplings 
of the Wess-Zumino-Witten Lagrangian [42, will also contribute to the spin polarisabilities. These terms are 
described below. 
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B. Baryons 



In SU(4|2) HBxPT, the physical nucleons (those composed of three valence quarks) are embedded in a 70- 
dimensional representation of the flavour group described by a three index flavour-tensor, B [23, US- Since the 
mass-splitting between the nucleon and A-isobar, A = Ma — Mjy ~ 300 MeV, is comparable to the physical pion 
mass (and less than pion masses used in current lattice simulations), the A-isobar must be included in the the- 
ory. These fields are represented in a totally symmetric three index flavour-tensor T M (a Rarita-Schwinger field) 
transforming as a 44-dimensional representation of SU(4|2). The mass-splitting A is small compared to the chiral 
symm etry breaking scale, and in this work we treat A ~ m, in the power counting [30j . For additional details, see 

Refs. Hamm. 

The relevant part of leading-order Lagrangian describing these baryons and their interactions with Goldstone mesons 

is 



C ( z ] = t(Bv VB) - i (T"v ■ VT^j + A (t% 



+2a (BS^BA^) + 2/3 (BS^A^B) + 2H (t* 'S» A^ + j|c \[t u A V B} + (BA U T)\ , (20) 

where is the baryon velocity, is the covariant spin-vector [2^, |^| and is the chiral, and U(l) covariant 
derivative M The axial-vector current is given by 

A» = l - [t;(d>*-ieQA>*)tf - g {d* -ieQA»)t] . (21) 

The various flavour contractions (indicated by the parentheses) are defined as in Ref. (3g . In order for T M to correctly 
describe the spin-3/2 sector, the constraints v ■ T = S ■ T = are used. The partially quenched Lagrangian, Eq. i|20[) . 
contains one more operator than the corresponding two flavor %PT Lagrangian. To determine the relation of the 
partially quenched operators to those of SU(2) x?T, one simply restricts the flavor indices of all the operators in 
Eq. (|2()|> to the valence sector, leading to 

4 1 Q 2 1 

C = (Jan , H = 3aa- (22) 

Here, g± is the coupling of the nucleons to the SU(2) singlet-meson field, which decouples in xPT. Thus when the 
QCD limit of the partially quenched theory is taken, where nij = m u and mi = mj, all dependence on g\ vanishes 
(in the quenched theory, quantities may still depend on this parameter, see Appendix lA"|l . 

As with the mesons, at leading order the photon is minimally coupled to the baryons with fixed coefficients. At the 
next order in the expansion there are a number of new electromagnetic gauge invariant operators which contribute to 
the Compton amplitude and the polarisabilities. Here, we display the relevant terms at this order, 

4 1} = H^ F ^ W (P P". S "\ + H P", S v \ Q i+ B) + /x 7 str [Q €+ ] (B [S», S»] B)} 

2^ T 2M^ [(^2«+ T '0 + (T'S v Q i+ B)] , (23) 

whe re U n.ri.-y are magnetic moment coefficients [3fil l38j. fj>x is the coefficient of the Ml transition 70-44 operator 
[HlH and 

Q s± - \@QZ±£Q&). (24) 

The partially quenched magnetic moment coefficients are related to the isoscalar and isovector magnetic coefficients, 
Ho and //i, in standard two flavour xPT as 



Mo 



+ M/3 + 2/i 7 ), Mi = g( 2 Ma ~ Vp)* ( 25 ) 

where the xPT Lagrangian describing the magnetic moments of the nucleons [the proton and neutron magnetic 
moments are /Lt Pi71 = ^(mo ± is given by 

F lw (iM,N[Si',S v ]N + iiJf[Si',Sr]>T% + N), (26) 



2M N 
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FIG. 2: Anomalous contributions to the polarisabilities. The crossed circle corresponds to the insertion of an operator from 
the Wess-Zumino-Witten Lagrangian, Eq. 1321 . and the crossed meson line corresponds to a hairpin interaction 



forr j± = 5 (£ + T°£±£"r a £ + )- 

There are other operators formally at this order which do not contribute to the polarisabilities at the order to which 
we work. There are kinetic operators and higher dimensional couplings of the baryons to the axial current whose 
coefficients are exactly fixed by the reparameterisation invariance of the baryon four-momentum |45l l4r| . These 
operators give the Z dependent pieces of the Compton amplitudes in Eq. J2J). There are also additional operators 
with unconstrained coefficients such as (B A - AB) that contribute to the Compton amplitude at higher order. In 
two-flavor x?T there are two such operators, and in the SU(4|2) case there are ten |46[ . 

The leading operators which contribute to the electromagnetic polarisabilities at tree level occur at 0(Q A ) and are 
given by the general form, 

p 2 p p p 

(gr""Q| ± g), 

x 

(where the T^ v are spin structures) while the leading tree-level contributions to the spin polarisabilities occur at 
(D(Q 5 ). The complete set of such operators in the case of two-flavour x?T is given in Ref. |JjJ. Again there are 
significantly more such operators in QxPT and PQxPT. We do not explicitly show these operators, as they do not 
contribute at the order we are working and will not modify volume dependence until 0(Q e ). 



V. NUCLEON POLARISABILITIES 



Using the Lagrangian of the preceding section, we can calculate the amplitudes defined in Eq. for Compton scat- 
tering from a nucleon (extensions to full octet and decuplet of baryons are straight-forward although the convergence 
of HB^PT with three-flavours is not clear) . We work with a power counting such that 

Q~e~|!~^~f (27) 
A x A x A x 

(it is also convenient to count A/A x as the same as Q as it is numerically similar at the masses relevant for current 
lattice calculations). Below, we will also restrict ourselves to the low frequency limit ui <C in order to extract 
the polarisabilities from the Compton scattering amplitudes defined in Eqs. Q and For larger energies, the 
concept of polarisabilities breaks down and the target essentially becomes a dispersive medium. Working to order 
Q 3 in the chiral expansion, Compton scattering requires the calculation of the diagrams shown in Figs. |2 [3] and 0] 
(and a corresponding set involving internal 44-plet baryons). By definition, tree level contributions from nucleon pole 
diagrams do not contribute to the polarisabilities; their contribution to the amplitudes are given explicitly in Eq. J2J. 
For each polarisability X = a, j3, 71-4, it is convenient to separate the different contributions as 

x = ^anomaly + X A + X l °° P , (28) 

corresponding to the contributions from Figs [21 121 an d El respectively. We discuss these contributions in the following 
subsections. At order Q 3 , all contributions are expressible in terms of a small set of LECs that contribute in many 



4 Loop and pole |4Sl contributions with 44-plet intermediate states must be included since A is a small-scale. Any A dependent terms 
analytic in ra^ arising from the loop diagrams, and additional operators pr oportio nal to powers of A/A x can be resummed into the 
appropriate LECs of A independent operators (the LECs then depend on A) I4(ll4ll . Keeping these contributions explicit is redundant 
as A can not be varied in a controlled manner. 
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FIG. 3: Born diagrams involving internal 44-plet states that give contribution to the polarisabilities. 




FIG. 4: Diagrams contributing to nucleon polarisabilities at order Q 3 . The solid and dashed lines denote nucleons, and pions 
respectively. Additionally, our results include a similar set of diagrams in which the internal 70-plct propagator is replaced by 
a 44-plet resonance. 

other processes and are thus reasonably well determined (at least in the xPT case). The total 0(Q 3 ) loop contribution 
is finite, but loop-contributions at higher orders are divergent; as discussed in the preceding section, the counter-terms 
specific to Compton scattering that absorb these divergences and the associated scale dependence enter at 0(Q 4 ) for 
the electric and magnetic polarisabilities and 0(Q 5 ) for the spin polarisabilities. 



A. Volume independent contributions to polarisabilities 



1. Anomalous contribution to "/N — > jN: ir — > 77 



The anomalous decay of flavour neutral mesons to two photons |5lL l52| has important consequences in Compton 
scattering in non- forward directions. These contributions arise from the meson pole diagram shown in Fig. EI Anoma- 
lous decays are well understood in x?T, entering through the Wess-Zumino-Witten (WZW) Lagrangian p2l |43| . 
However, these effects have not been investigated in the quenched and partially quenched theories and some interest- 
ing subtleties arise. 

In SU(2) xPT, the one-pion, two-photon piece of the WZW Lagrangian, is given by 



where 



3e 2 



16tt 2 / 



rtr 



uvpo p T? 



[iv -l pa 1 



(29) 




(30) 
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This Lagrangian is completely determined as its coefficient can be fixed by directly matching to the perturbative 
QCD calculation of the relevant triangle diagram (the one loop calculation is exact [55j, in accordance with Witten's 
geometric quantisation condition |43|). At higher orders, additional anomalous operators appear [3] but they do not 
contribute to Compton scattering until 0(Q 5 ). 

It is well known that quenched and partially-quenched chiral perturbation theories generally have more complicated 
operator structure than in the case of QCD (e.g., one can not use Cayley-Hamilton identities 55]). Thus, in order to 
generalise Eq. I|29|> to the partially-quenched cases, we might imagine the extended tt° — > 77 Lagrangian to be of the 
form 



r p Q nr f f"^p<? p p 

7r°77 V P a 



ax str [$Q 2 ] + a 2 str [$Q] str [Q] + a 3 str [$] str [Q] 2 + a 4 str [$] str [Q 2 



(31) 



(in the quenched case only the first operator is non- vanishing, but a similar discussion applies). With the condition that 
in the QCD limit where the sea-quark and ghost-quark masses and charges are set equal to those of the valence-quarks, 
matrix elements of Eq. I|31ll reproduce the matrix elements of Eq. I|29() . As discussed in Section llVl the singlet field, $0 
acquires a large mass from the strong U(1),a anomaly 39] and can be integrated out of the partially-quenched theory; 
consequently, the operators proportional to 03 and 04 can be ignored. Additionally, from the multiple super-trace 
structure, one can deduce that the operators 

str[$Q]str[Q] , str [$] str [Q] 2 , str [$] str [Q 2 ] , 



have at least two closed loops at the quark level. Following the arguments in Refs. |51ll52ll53| . one can show that these 
operators do not correspond to anomalous quark level processes. Moreover, the leading dependence of the underlying 
quark-level diagrams is proportional to the quark mass, and thus the coefficient of these operators must scale as, 
^2,3,4 ~ Tiq/A 2 ,. Although these operators contribute to r\ a — ► 77, they are not anomalous, and only contribute at 
higher orders in the chiral expansion. We can thus conclude that the only operator in the Lagrangian, Eq. 131f) . 
which contributes to the anomalous decay of the neutral mesons at leading order is str [<I>Q 2 ] . The coefficient is easily 
determined by matching to either perturbative partially-quenched QCD or to the xPT expression in the QCD limit. 5 
The appropriate Lagrangian is therefore 



3e 2 



16tt 2 / 



str [$Q 



21 ^flVpCT 



F F 



pa 



(32) 



From the above Lagrangian, it is apparent that all of the flavour diagonal fields in Eq. (|14|1 . have anomalous couplings 
to two photons. Calculating the diagrams in Fig[2|leads to the following anomalous contribution to Compton scattering 
on the proton 6 in partially- quenched 



rpPQ, anomaly 



i tiivaft k kP v • S 



91 



24e 2 I 



1 2 



UU 

'2 



2 qi 

1 2 

2* 



1 



ll. 



2 qi 



(1] flf ) 



'■(iii 



(r 2 - m\ u )[r 2 - m\) 



(r 2 — m 2 x ) \ (r 2 — m 2 M ) (r 2 



"ddi 



(33) 



In the above expression, r = q 1 — q, is the momentum transfer to the nucleon and A^- = m 2 z — m 2 ^- is a measure of 
the isospin breaking in the sea sector. In the sea isospin limit (mi — > m 3 ), the double pole structure of the amplitude 
vanishes, and in the QCD limit all dependence on g\ vanishes. 

Expanding Eq. 133fl in frequency and comparing with Eq. Q leads to the following anomalous contributions to the 



5 We can thereby bypass the need to extend Witten's global quantisation condition |4.'I to non-compact graded Lie groups. 

6 The anomalous contribution to neutron-Compton scattering is given by Eq. 1331 with the interchange of u <-> d. 
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polarisabilities: 



^anomaly _ 
□anomaly 



anomaly 

7i 

anomaly 
72 

anomaly 
7 3 



anomaly 
74 



3e 2 G a 



87T 3 / 2 77l2 



= 0. 



3e 2 G a 



167T 3 f 2 m 2 ' 
3p 2 G 



(34) 
(35) 

(36) 

(37) 

(38) 

(39) 



167T 3 / 2 TO 2 ' 

where the coefficients, G an om, are given in Tabled for the different theories under consideration. These contributions 
vanish in the iso-scalar combination of proton and neutron targets in the QCD limit. 



2. A resonance contributions 



The contributions to the amplitude from the Born-terms involving the 44-plet resonance (which contains the A- 
isobar), Fig. EI are identical in xPT, PQxPT and QxPT as they are purely valence quark processes. They are given 
by 



£ 



7i 
7 2 A 



7 3 



74 



= 0. 



/4 



e 2 (q u 



367r(2AfAr) 2 A ' 



0, 



2 e {Qu 



Id) 



72tt(2M n ) 2 A 2 



„2 e2 (lu 



Qd) 



72n(2M N ) 2 A 2 ' 

where /j,t is the magnetic dipole transition coupling of Eq. i|23|) . 



(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



B. Infinite volume 



The loop contributions to infinite volume chiral expansion of the polarisabilities in %PT are well known at order 

o 3 ^iliiiilelilil and at o(q 4 ) mmililililiilz^ ^ at this order ' the A_resonances have not 

been included as dynamical degrees of freedom, restricting the range of applicability to m v <C A). Since the photon 
only couples to charged mesons, the results in the quenched and partially-quenched theories are similar to those 
in xPT. In particular, no quenched or partially-quenched sicknesses (double pole contributions from neutral meson 
propagators) enter expression for the loop diagrams. In general, the quenched power counting presents differences for 
electromagnetic observables [33l |36l 13 UM > however no new contributions appear at the order we work. 

Using the effective couplings Gb , G' B , Gt and G' T given in Table [IJ we find that the loop contributions to the 
polarisabilities are 



loop 



^loop = 



e 2 


'5G B 1 


5G' 1 
h B 


Gt 


4^/ 2 


1927T 


1927T m u j 


F 727T 2 


e 2 


' G B 1 


G' 1 

V B 


Gt 


4^/ 2 


384-7T 


3847T m u j 


h 72tt 2 



F Q (m w , A) 



72^r 



G' 



A) 



Fp(m n ,A) + 7 ^F f3 (m uj ,A) 



(46) 



(47) 



7i 



loop 



4tt/ 2 



G B 1 
487r 2 ml 



+ 



G'r 



48^ 2 m% 



Gt 
72^ 



72tt 



-F 71 (m Mj , A) 



(48) 



13 



QCD 



QQCD 



PQQCD 



gA^Zql + {1 - Z)q 2 d ^ - g? - ft 2 ] 
+5i {id ~ <7j - if + <lu) 



G anom g A (2Z - l)(g 2 - g 2 ) 2g A (Zq 2 u + (1 - Z)q 2 d ) 



+5i (3d + ql) 



G B 
G' B 
Gt 

Gyp 



4g 2 A (qd 



tffAiV (id - qu) 



§ (4ffi - 45A3i - 5# 2 ) (g d - g u ) 2 







-| (5p 2 + 4c/a5i - 4gi) (g d - g u ) 2 

| [ (6<?1 - 6 + qi) q d + 5g 2 + 5qf + 4g 2 - 4g j g„ - 4g ; g„) # 2 
+4ffA (g 2 - 2g u g J - + g 2 + 2ql - 2g ; g„) ff i 
+8ffi (g 2 - 2g„g 3 + qf + 2g 2 - 2g ; g„) ] 

5 2/ \2 

s5an Wd - qu) 

San (4gJ - 4 (g 3 + g ; ) g d + 3g 2 + 3gf + 2g 2 - 2q 3 q u ~ 2q t q u ) 



1 „2 
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TABLE I: Effective couplings for the various contributions to the polarisabilities. 



loop 
72 



47T/2 



G" R 



967r 2 to 2 967T 2 ra 2 - 



2 ^72 

TYl u j , 

A) 



72tt 



72tt 



(49) 



loop 

7 3 = 



4^/2 



G B 1 
192tt 2 m 2 



Go 



192tt 2 m 2 



^^(m^A) + Jg^^, A) 



(50) 



where 



loop 

74 



F a (m,A) 
F (m,A) 
F 71 (m,A) 
F 72 (m,A) 



67 



G' 



192vr 2 m 2 192tt 2 m 2 ■ 



9A 



10m 2 



2(A 2 - to 2 ) 3 / 2 



In 



Va 2 



rrir + «e 



VA 2 



m' + ze 



1 



2(A 2 -m 2 ) 1 / 2 
A 2 + 2m 2 



In 



A - VA 2 - m 2 + ie 
A + a/A 2 " 

3Am 2 



m + ie 



(A 2 -m 2 ) 2 2(A 2 -m 2 ) 5 / 2 
F 73 (m,A) =F 74 (m,A) 



In 



A - VA 2 - m 2 + ie 



VA 2 



1 



m 2 + ie 



In 



m 2 2(A 2 -m 2 ) 3 / 2 



A - \/A 2_ 



m 2 + ie 



A + VA 2 - to 2 + ie 



(51) 

(52) 
(53) 
(54) 
(55) 



Here we have used dimensional regularisation, however the results are finite and hence independent of the regulator 
without the addition of counterterms. These loop contributions vanish at zero photon frequency, preserving the 
Thompson limit. They are identical for both proton and neutron targets, though isospin breaking effects from loops 
enter at 0(Q 4 ) in the expansion. In the x?T case, these results reproduce those of Refs. |5^.lr30| . 



Finite volume 



In momentum space, the finite volume of a lattice simulation restricts the available momentum modes and conse- 
quently the results differ from their infinite volume values. These long-distance effects can be accounted for in the 
low-energy effective theory. Here we shall consider a hyper-cubic box of dimensions L 3 x T with T L. Imposing 
periodic boundary conditions on mesonic fields leads to quantised momenta k = (fco, k), k = =j£-j = ^-(ji, J2, J3) with 
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ji G Z, but ko treated as continuous. On such a finite volume, spatial momentum integrals are replaced by sums over 
the available momentum modes. This leads to modifications of the infinite volume results presented in the previous 
section; the various functions arising from loop integrals are replaced by their finite volume (FV) counterparts. In a 
system where m^L 3> 1, the power counting of the infinite volume low-energy effective theory remains valid and finite 
volume effects are predominantly from Goldstone mesons propagating to large distances where they are sensitive to 
boundary conditions and can even "wrap around the world" . Smaller volumes in which m^L ~ 1 are discussed in 
Appendix IbI Since the lowest momentum mode of the Goldstone propagator is ~ exp(— m^L) in position space, finite 
volume effects will behave as a polynomial in 1/L times this exponential if no multi-particle thresholds are reached 
in the loop (as is the case in these calculations provided the photon energy is small enough, u < m v ). 

Repeating the calculation of the loop diagrams using finite volume sums rather than integrals leads to the following 
expressions for the loop contributions to the polarisabilitics: 



loop 



(L) 



1152tt/ 2 



dX 



ZG B T a {M uu ) + 3G' B f a (M uj ) + 8G r ^(M„ B ) + 8G' T F a (M*) 



(56) 



(3 loop (L) 



1152tt/ 2 



dX 



SGbMMuu) + 3G'MM uj ) + 8G T Fp(Mt) + 8G' T F P (M%) 



(57) 



7i oop (£) 



7e 2 



576tt/ 2 



dX 



3G B ^ 71 (M uu ) + 3G^ 7l (M uj ) - 4G T ^" 71 (M* ) - 4G^ 7l (M 



(58) 



64^/ 2 jo 



dX 



3G B F l2 (Muu) + 3G' B F l2 (M uj ) - AG T T l2 {Mt u ) ~ ±G' T T l2 {M%) 



(59) 



7e 2 

loop / j \ 1 ° 

73 ^ ~ 11527T/2 



3G B F l3 {M uu ) + 3G' B ^ 3 {M uj ) - 4G T ^ 3 {M^ U ) ~ ±G' T T l3 (M 



(60) 



7 r p w = wr p (L), 

where M ab = \Jm\ h + A 2 and = yjm 2 ab + 2AA + A 2 and 



(61) 



T a (m) 


= 180A 2 Tt(to 

2 v 


) + 190 Jt (to 


- 280A 2 Jo (m) - 455/C 9 (to) + 315A 2 /C ii (to) + 252£ ii (m) , 


(62) 


Tp(m) 


= 60,77 (rn) - 

2 v ' 


224/C 9 (to) + 

2 v ' 


189£h(to) , 

2 v ' 


(63) 


Fyi M 


= 30A 3 X 9 (to) 

2 v ' 


+ 10A J 9 (to) 

2 v ' 


- 45A 3 Jii (to) - I8A/C41 (to) , 


(64) 


T l2 (to) 


= XJC n (to) , 






(65) 


T 13 (to) 


= 10 A^ (to) - 

2 v ' 


- 9A/C ii (to) , 




(66) 



and the finite volume sums Xg(m), . . . , Cp{m) are defined in Appendix IH1 These expressions reduce to the results of 
Eqs. (|46[) - l|51[) above in the infinite volume limit. 

To illustrate these effects, Figs. [5] and HO show the volume dependence of the various polarisabilities in the proton 
and the neutron, respectively. Here we have specialised to QCD, setting q u = 2/3, qd = —1/3, §a — 1-25, |<?ata| = 1.5, 
fix = 10.9, / = 0.132 GeV, M N = 0.938 GeV and A = 0.3 GeV. 7 In each plot we show results for the ratio 



AX(L) 



X(L)-X(oo) 
Xjpo) 



(67) 



for the six polarisabilities at three different pion masses, m n = 0.25, 0.35, 0.50 GeV. The overall magnitude of these 
shifts varies considerably; generally volume effects are at the level of 5-10% for = 0.25 GeV and smaller for larger 



7 The value of (it is chosen to correspond to that found in analysis of Ref. |6fl {fix = 2\/2b\ of that reference). In principle this LEC 
can be determined from an analysis of lattice polarisabilities or TV- A transition matrix elements. 
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L [fm] l [fm] 



FIG. 5: Volume dependence of the proton polarisabilities. Here we show the ratio of the difference of the finite and infinite 
volume results to the infinite volume results for three values of the pion mass using the parameters described in the text. The 
curves terminate at m„ L = 3. 

masses. Larger effects are seen in a number of the spin polarisabilities but there are as yet no lattice calculations of 
these quantities. The magnetic polarisability has a particularly small volume dependence which can be understood 
from the large decuplet resonance contribution that is independent of the volume. 

The above expressions also allow us to calculate the finite volume effects in the quenched data on the various 
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m^=250 MeV 
m T =350 MeV 
m n =500 MeV 
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FIG. 6: Volume dependence of the neutron polarisabilities. The various curves are as in Fig. 



polarisabilities calculated in Refs. 0,0- The quenched expressions involve a number of undetermined LECs (quenched 
9 A, <?ij 9na and /it are unrelated to their PQ^PT/xPT values), so we can only estimate the volume effects. To do 
so, we choose q u = 2/3, q d = -1/3, g A = 1.25, gi = 1, \g NA \ = 1.5, \i T = 5.85, / = 0.132 GeV, M N = 0.938 GeV and 
A = 0.3 GeV and take the pion masses corresponding to the lightest used in these lattice calculations, ~ 0.5 GeV 
(we ignore the issue of the convergence of xPT at such masses) . The results for the volume dependence of the various 
polarisabilities of the proton and neutron are shown in Figs. and [SJ In each plot, the shaded region corresponds 
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FIG. 7: Volume dependence of the proton polarisabilities in quenched QCD at the lightest quark mass used in the lattice 
calculations of Refs. 0,13. The central curves and shaded region correspond to the parameters quoted in the text. 



to reasonable variation of the unknown couplings, — 1 < g\ < 1, 0.8 < |.9ata| < 2 and 2.8 < \ht\ < 8.5. From these 
figures, we see that the calculations on a (2.4 fm) 3 lattice with =0.5 GeV may differ from their infinite volume 
values by 5-10% in the case of the electric polarisability and a few percent for the magnetic and spin polarisabilities. 

An interesting effect that arises at finite volume is that the Thompson-limit and other Born terms in the frequency 
expansion of the scattering amplitude (terms in Eq. (|5J that are not polarisabilities) receive finite volume contributions 



18 





0.002 




£ -0.002 



-0.004 
-0.006 



4 

L[fm] 






006 





004 





002 







-0 


002 


-0 


004 


-0 


006 



4 

L[fm] 




FIG. 8: As in Fig. but for the neutron. 



from the loop diagrams in Fig. 0] that vanish exponentially as the volume is increased. As an example, the amplitude 
for Thomson-limit (zero frequency) scattering on the neutron (which is identically zero at infinite volume) is shown 
in FigEI This result is somewhat counter-intuitive, but arises from the effects of the periodic boundary conditions on 
the long range charge distribution of the hadron. It does not imply the non-conservation of charge. 

The results presented here all assume that the higher order terms in the Q expansion provide small contributions 
to the volume dependence of the polarisabilities. This may or may not be the case as diagrams that are formally 
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L[fm] 

FIG. 9: Volume dependence of the Thompson limit of photon neutron scattering. Notice the infinite volume limit is zero. 

of higher order in the infinite volume vPT power-counting can have volume effects that are enhanced over those at 
lower infinite volume order (see Ref. 41] for a detailed discussion). Such issues may be particularly relevant for the 
polarisabilitics where the convergence of the chiral expansion is tenuous. In this regard, studying the FV behaviour of 
the lattice results may in fact be a useful diagnostic tool with which to determine if or why the convergence is poor. 

VI. CONCLUSION 

We have investigated Compton scattering from spin- half targets from the point of view of lattice QCD. We first 
discussed how external field methods can be used to probe all six polarisabilities of real Compton scattering for both 
charged and uncharged targets. Such calculations will tell us a lot about the low energy QCD structure of hadrons 
and will be of much use in phenomenological studies requiring the full set of polarisabilities as only certain linear 
combinations are available from current experiments Q, . The techniques discussed here also allow us to extract 
other electric properties of charged particles using external fields including the electric dipole moment of the proton 
and the quadrupole moment of the deuteron. 

Our second major focus was on the effects of the finite volume used in lattice calculations on the polarisabilitics. 
Since polarisabilities are infrared-sensitive observables (they scale as inverse powers of the pion mass near the chiral 
limit), the are expected to have strong volume dependence. This is indeed borne out in the explicit calculations 
presented here. In QCD, we generically find that the polarisabilities experience volume shifts of 5-10% from the infinite 
volume values for lattice volumes ~(2.4 fm) 3 and pions of mass 0.25 GeV. The electric and first spin polarisabilitics 
are particularly sensitive. In the case of quenched QCD (relevant to the only existing lattice data), we find significant 
effects even at pion masses ~ 0.5 GeV. Future lattice studies of the polarisabilities should take these effects into 
account in order to present physically relevant results. 

As extensions of this work, one can also consider the generalised polarisabilities, higher-order polarisabilities and 
parity violating polarisabilities (see Ref. [T^) all of which can be extracted from appropriate lattice calculations 
similar to those detailed in Section 11111 Such information would lead to a further-improved understanding of the 
low-energy structure of the hadrons and prove very useful in directing the next generation of precision Compton 
scattering experiments. The lattice provides a novel opportunity to study the neutron polarisabilities directly instead 
of from nuclear targets and extending the lattice methods of Section ITTT1 to the deuteron (along similar lines to those 
discussed in Ref. |73|) will also prove useful for comparison to experiment. 
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APPENDIX A: QUENCHED CHIRAL LAGRANGIAN 

In this Appendix, we display the relevant pieces of the quenched chiral Lagrangian in the meson and baryon sectors 
and note particular pathologies of the quenched theory In a quenched two flavor theory, we have valence (u, d) and 
ghost (£t, d) quarks with masses contained in the matrix 

m Q = diag(m„, md,m fi ,m^) , (Al) 

where m~ j = m u _d to maintain the exact cancellation from the path-integral determinants arising from the valence 
and ghost quark sectors. The corresponding low-energy meson dynamics are described by the Q%PT Lagrangian. 
At leading order, the form of this Lagrangian is the same as in Eq. (|12fl where the pseudo-Goldstone mesons are 
embedded non-linearly in £ with the matrix $ now given by 



*= ?- > (A2) 





where 

M= [ nu _ *\ , M= I ^ * + ) , x=\ AVu A " T ) • 

\7T Vd J \7T fj d J 

The matrix M contains the usual valence-valence mesons, while mesons in M are composed of ghost quarks and 
anti-quarks, and finally those in \ of ghost-valence quark-anti-quark pairs. Unlike the partially quenched theory, 
there is no strong U(l),4 anomaly, and the flavor-singlet field, $o = str[$]/\/2 (along with its couplings mo and a$), 
must be retained in the theory. For the electromagnetic and spin polarisabilities in QQCD, no loop contributions 
from the singlet are needed to the order we work as flavor-neutral mesons are not present in loop diagrams at this 
order. Despite flavour neutral mesons being absent in loop graphs, the anomalous tree-level term couples the quenched 
singlet to the nucleoli. Cancellations, however, lead to final results that are independent of mo and a$. 
For the quenched electric charge matrix of the valence and ghost quarks, we choose 

Q = dimg(q u ,q d ,q u ,q d ) . (A4) 

Notice the peculiarity that str Q — is unavoidable in the quenched theory. In the quenched theory, there are 
anomalous decays of flavour neutral mesons into two photons. In terms of SU(2|2) QQCD quark fields, contributions 
to the anomaly from the valence and ghost sectors come weighted with squares of the quark charges, and we are thus 
not restricted to only the flavor singlet current (as is the case for the strong U(1)a anomaly). The relevant term of 
the anomalous quenched chiral Lagrangian is the same as has been detailed above in Sec. IV ATI 

In SU(2|2) HBxPT, the nucleons (those composed of three valence quarks) enter as part of a 20-dimensional 
representation described by a three index flavour-tensor, B. The quenched A-isobar is contained in the totally 
symmetric three index flavour-tensor transforming in the 12-dimensional representation of SU(2|2). The leading- 
order Lagrangian describing these baryons and their interactions with Goldstone mesons is 

Cf q = i(Bv VB) + 2a (~BS»BA„) + 2/3 (BS^A^B) + 2 7 (BS^B) stiv^ (A5) 
-i (T% ■ VT^ + A (rX) + 2H (T U S»A^ + 2 7 ' str^ + J^C [(T ' A v B) + (~BA„T» 

In contrast to partially-quenched and unquenched chiral perturbation theory, there are two additional axial couplings 
7 and 7' due to the presence of the flavour-singlet field. One should keep in mind that although we use the same 
notation for simplicity, all of the coefficients in the quenched Lagrangian have distinct numerical values from those of 
the partially-quenched Lagrangian. In the large N c limit, the coefficients of the two theories are related [74j . 
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Again the photon is minimally coupled in the above Lagrangian with fixed coefficients. At the next order in the 
expansion, the relevant terms that appear are 

+V|mt^^ [(BS^Qt+T") + (T^ S u Q.£ + B)] . (A6) 

The PQxPT term with coefficient /x 7 is absent in the quenched theory. This only affects the Born terms of the 
Compton amplitude, which are essentially unknown because they depend on the quenched magnetic moment. Finally, 
the leading two-photon operators that give completely local contributions to the Compton scattering tensor appear 
in quenched chiral perturbation theory in essentially the same form as PQxPT. However, there are fewer operators 
per spin structure compared to the partially quenched case because of the super-tracelessness of the electric charge 
matrix. Our computation is unchanged since these terms do not contribute at the order we work. 



APPENDIX B: COMPTON SCATTERING IN SMALL VOLUMES 



In this appendix, we discuss the volume dependence of the polarisabilitics on asymmetric lattices in which the 
spatial dimensions are small but the temporal direction remains large: m f L <C 1 but L 3 Tm q (qq} 3> 1. In this 
e'-regime [7(|, m w ~ e' 2 and L ~ e' _1 (where e' is the small expansion parameter) so Goldstone boson zero modes 
(modes with q = Q) are enhanced, but remain perturbative. Thus, the power counting of loop diagrams involving 
zero-modes is modified. In Compton scattering at finite volume, the photon momenta provide additional scales 
whose power-counting must be specified. The spatial components of the incoming and outgoing photon momenta 
are quantised as if- 1 ' — -£-fiqW (where ft q w are integer 3-tuples) and consequently scale as 0(e'). For real Compton 

scattering, qx^ 2 = |<7^| 2 , implying that the frequency u> is either zero or also 0(e'), parametrically larger than the 
pion mass, cj/m n 3> 1. This introduces an intrinsic difficulty in defining polarisabilities for realistic pion masses in 
small volumes, as they arise from an expansion around the zero frequency limit. Thus our discussion is restricted to 
the full amplitudes Ai(u>, 9) in Eq. QJ. This is not an issue in larger volumes fSec. IV 01 above ) as the quantisation of 
momenta is fine-grained on the scale of m^. In (doubly) virtual Compton scattering, the on-shell condition is relaxed 
and one can again consider non-zero energies comparable to or smaller than the pion mass. 

In the e'-regime, the loop diagrams in Fig. 0] generically contribute at order e' 3 for non-zero modes (we count 
e ~ e'). However, diagrams 3(a) an 3(b) contain non-derivative couplings and zero-modes provide a further enhanced 
contribution, ~ e'. Although the remaining diagrams in Fig. 0] contain derivatively coupled pions, the non-zero 
momentum insertions allow the energy integral to be performed with a pole ko ~ ~ e' 2 , leading to a putative 
enhancement. However, transversality of the Compton amplitude causes these enhanced contributions to vanish. 
Thus A±(u) = 0) receives contributions at O(e') and A3 at 0(e' 2 ) and these amplitudes will exhibit enhanced volume 
dependence for lattice calculations in the e'-regime. This dependence is given by 

A 1 ( U j 9) = e2 ° B 1 2e2Gr 1 1 <n{c' 3 ) (Bl) 

MUi6) = f GB " 2 -^%- ( ^ , + 0(e' 3 ). (B2) 

v ; 4m w L 3 m 2 -w 2 im^L 3 (m^ + A) 2 — lo 2 y ' y ' 

Since is necessarily smaller than any non-zero value of lo in the small volume regimes, pions can go on-shell in 
the loop diagrams of Fig. 21 Consequently the power-law dependence on volume seen in Eqs. idU and is not 
unexpected. The other amplitudes have contributions from the diagrams in Fig. 0] starting at C(e' 3 ). At this order 
additional diagrams such as those arising from tadpole dressings of the vertices in Fig.0Ja,b) also contribute and the 
full results arc left to future work. 



The calculation of the volume dependence of the polarisabilities in the e-regime |7.^ (where L ~ m^T ~ 1 and L 3 Tm q {qq) ~ 1 and 
zero-modes become non-perturbative) is beyond the scope of this work. 
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APPENDIX C: FINITE VOLUME FUNCTIONS 



The sums required in the evaluation of the polarisabilities at finite volume are (k = , with n a triplet of integers) 



MM) = 



g [|fc| 2 +M 2 ] 

Jfs(M) = Ip- X {M) - M%(M) , 

K,p(M) = Jg_ 2 (M) — 2M 2 J l g_i(M) + M A lp(M) , 

Cf)(M) = l^ z {M)-2,M 2 l fj ^ 2 {M) + 2,M A l^ 1 (M)- M e l fj (M) . 

At infinite volume these can be simplified using, 

i r(/?-§) 1 



lp{M,L^ oo) 



(4tt)§ r(/3) (M 2 )' 3 " 



for /? > 3/2. 

In numerically evaluating these sums, it is useful to note that 



£■ 



i 



^ (|n| 2 + x 2 )' 3 

n 

he expon 
exponentially fast in |n|. 



£ 



^^(Inp+x 2 ) 7T - 



r(/3) 



r(/3) Jo 



£ 



e + 1 



(CI) 

(C2) 
(C3) 
(C4) 



(C5) 



(C6) 



where E n (x) is the exponential integral function. This form is valid for [3 > |, x £ R and the remaining sums converge 
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